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y—\ , Abstract. We have collected the parameter of 38 neutron stars (NSs) in binary systems with 
spin periods and measured masses. By adopting the Boot-strap method, we reproduced the 
procedure of mass calculated for each system separately, to determine the truly mass distribution 
of the NS that obtained from observation. We also applied the Monte-Carlo simulation and 
introduce the characteristic spin period 20 ms, in order to distinguish between millisecond pulsars 
(MSPs) and less recycled pulsars. The mass distributions of MSPs and the less recycled pulsars 
could be fitted by a Gaussian function as 1.45±O.42M0 and 1.31±O.17M0 (with la) respectively. 
As such, the MSP masses are heavier than those in less recycled systems by factor of ~ 0.13Mq, 
since the accretion effect during the recycling process. 
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The measured mass distribution of NSs in pulsar binary systems is more diverse than 
previously thought. We have applied the Boot-strap method together with the Monte- 
Carlo simulation to the eclipse X-ray binaries, low mass X-ray binaries, double neutron 
stars (DNS) and white dwarf (WD)-NS systems, respectively. The result of this simulation 
t-H I is shown in Fig. [TJ Fitting a Gaussian function gives us the mass distribution is 1.31 ± 
0.29M Q with lcr confidence level (dash line). While the bimodal distributions can be 
fitted the simulated result perfectly. This result agrees with some theoretical predictions 
such as those by Kiziltan et al. (2010) and Valentim et al. (2011). The solid line shows 
the double Gaussian fitting curves (1.30±0.15M Q and 1.41±0.59M Q with la confidence 
level, respectively). It should be noticed that the bimodal distribution is coincidence 
with the result of Bayesian statistical method based on the data of DNS and NS-WD 
systems as showed by Kiziltan et al. (2010). This has strongly supported the idea that 
the recycled MSP become more massive (~ 0.13M Q ) than those in less recycled pulsar 
via the accretion phase in recycling process. To discriminate clearly whether this result 
is due to the selection effect or different evolution scenario, we set a characteristic spin 
period (P S pm = 20 ms) which helps us to divide all the sources into two groups: MSPs 
with P sp in < 20 ms and less recycled pulsars with P sp i n > 20 ms (zhang et al. 2011). 
We find that the average distribution of MSPs is 1.48 ± 0.53M Q while in less recycled 
pulsar is 1.35 ± 0.27M Q . The results are shown in Fig. [2] which indicate that the mass 
of rapidly rotating MSPs is systematic higher than those in the less recycled ones by 
factor of ~ 0.13Mq during the accretion phase. Furthermore, the large deviation of the 
MSP distribution could be reduced if we invoke the accretion-induced collapse (AIC) in 
WDs, since the AIC is expected to produce normal NSs, which in binaries can evolve 
into MSPs through the usual recycling scenario. 



1 



2 



Z. Chengg et al. 



0.7 



0.6 - 




NS mass[Msun] 

Figure 1. Monte Carlo simulation of NS masses 
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Figure 2. Monte Carlo simulation of mass for MSPs (Left) and less recycled pulsars (Right). 
2. Conclusion 

The simulation results of the average masses show that the MSPs are heavier than 
those in less recycled pulsars. This result strongly support the idea that the MSPs ac- 
creted ~ O.13M during the accretion phase in the recycling process. Furthermore, the 
AIC scenario must be invoked in order to reduce the large deviations of the MSP dis- 
tribution. This work is supported by NBRPC(2009CB824800, 2012CB821800) and the 
NSFC(10773017, 11173034). 
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